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Abstract. The method of solving the Bethe-Salpeter equation in Minkowski space, developed previously for 
spinless particles [T], is extended to a system of two fermions. The method is based on the Nakanishi integral 
representation of the amplitude and on projecting the equation on the light-front plane. The singularities 
in the projected two-fermion kernel are regularized without modifying the original BS amplitudes. The 
numerical solutions for the J = bound state with the scalar, pseudoscalar and massless vector exchange 
kernels are found. The stability of the scalar and positronium states without vertex form factor is discussed. 
Binding energies are in close agreement with the Euclidean results. Corresponding amplitudes in Minkowski 
space are obtained. 
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1 Introduction 

Bcthc-Salpctcr (BS) equation for a rclativistic bound sys- 
tem was initially formulated in the Minkowski space [2]. 

' It determines the binding energy and the BS amplitude. 
However, in practice, finding the solution in Minkowski 
space is made difficult due its singular behaviour. The sin- 
gularities are integrable, but the standard approaches for 

• solving integral equation fail. To circumvent this problem 

' the BS equation is usually transformed, by means of the 
Wick rotation, into Euclidean momentum space. 

The binding energy provided by the Euclidean BS equa- 
tion is the same than the Minkowski one. For computing 

■ the binding energy it is enough to solve the Euclidean BS 
equation in the rest frame p = 0. While the rest frame 
solution is not enough to obtain the electromagnetic form 
factors, since the initial and final states cannot be simul- 
taneously at rest. Contrary to the Minkowski amplitude, 
the dependence on the total momentum p of the Euclidean 
amplitude is not given by a standard boost, but must be 
found numerically by solving the Euclidean BS equation 
for non-zero p. This requires the rather complicated cal- 
culations which were performed in [210]. Within this ap- 
proach, the hadronic form factors have been calculated, 
albeit with the aid of additional assumptions. 

However, the knowledge of the Euclidean amplitude in 
a moving system is still not sufficient to calculate some ob- 
servables, e.g. electromagnetic form factors. The integral 
providing the form factors contains singularities which are 
different from those appearing in the BS equation. Their 
existence invalidates the Wick rotation, to the zero com- 



ponent ko , in the form factor integral [3] and prevents from 
obtaining the exact form factors in terms of the Euclidean 
BS amplitude alone. To avoid this problem, the knowledge 
of the BS amplitude in Minkowski space is mandatory. 
Thus, fifty years after its formulation, finding the BS so- 
lutions in the Minkowski space is still a field of active 
research. These solutions would pave the way to numer- 
ous applications going from the hadronic electromagnetic 
form factors of mesons to the deuteron electrodisintegra- 
tion amplitudes. 

Some attempts have been recently made to obtain the 
Minkowski BS amplitudes. The approach proposed in [S] 
is based on the integral representation of the amplitudes 
and solutions have been obtained for the ladder scalar case 
[51171IH] as well as, under some simplifying ansatz, for the 
fcrmionic one [5]. Another approach [TU] relies on a sepa- 
rable approximation of the kernel which leads to analytic 
solutions. Recent applications to the np system can be 
found in [TT] . 

In a previous work [1] we have proposed a new method 
to find the BS amplitude in Minkowski space and applied 
it to the system of two spinless particles. 

Our approach consists of two steps. In the first one, 
the BS amplitude is expressed via the Nakanishi integral 
representation [T21I13] : 



$(k,p) = dz' dry' (1) 
J -X Jo 

v gJM 

X q" ■ 

[k 2 + p ■ k z' + \M 2 - m 2 - 7' + ie] ' 
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Notice that in this representation, the dependence on 
the two scalar arguments ft 2 and p-k of the BS amplitude 
is made explicit by the integrand denominator and that 
the weight Nakanishi function (7(7, z) is non-singular . By 
inserting the amplitude (p} into the BS equation on finds 
an integral equation, still singular, for g. 

In the second step, we apply to both sides of BS equa- 
tion an integral transform - light-front projection [T] - 
which eliminates singularities of the BS amplitude. It con- 
sists in replacing k — > k + ^ (3 where lu is a light-cone 

four- vector u> 2 = 0, and integrating over /3 in infinite lim- 
its. We obtain in this way, an equation for the non-singular 
g("f, z). After solving it and substituting the solution in eq. 
JTJ), the BS amplitude in Minkowski space can be easily 
computed. 

As a first application we have considered the spin- 
less case with ladder kernel [1] . The binding energies were 
compared with the direct solution in the Euclidean space 
and found to agree each other with high accuracy. By in- 
serting the computed weight function g in ([T]) and set- 
ting fco = ik^ the result was found to coincide with the 
corresponding Euclidean BS amplitude, found indepen- 
dently. The method has been also successfully applied to 
the cross-ladder kernel [T3] and to compute the electro- 
magnetic form factors [5]. 

The main difference between our approach and those 
followed in [6l[Tj[8l[9] is the use of the light-front projection. 
This eliminates the singularities related to the BS Minko- 
wski amplitudes. The method is valid for any kernel given 
by the irreducible Feynman graphs. 

We present in this paper the extension of our preced- 
ing work [1] to the two fermion system. In this case the 
Nakanishi function g is replaced by a set of n functions 
gi, satisfying a system of coupled integral equations, with 
n depending on the total angular momentum of the state 
(e.g. n~A for J = 0, n=8 for J = 1). We will see that the 
direct application to the fermionic kernels of the method 
used in the spinless case, is however marred with some 
numerical difficulties. Although they could be overcome 
by a proper treatment of the singularities, in this work we 
propose an alternative method allowing to solve the BS 
equation for two fermions in Minkowski space with the 
same degree of accuracy than for the scalar case. The nu- 
merical applications will be limited to the J 71 " = + state. 

The system of equations for the Nakanishi weight func- 
tions gi is derived in sect. [2] starting from the original 
BS equation. In sect. [3] we develop a regularization proce- 
dure for fermionic kernels. Numerical results for the scalar, 
pscudoscalar and massless vector exchange couplings are 
presented in sect. |H Section \b\ contains concluding re- 
marks. Some details of the calculations are given in the 
appendices [21 [B] and [Cj 



2 Derivation of equation 

The BS equation for the two fermions vertex function r 
reads: 

/d 4 k' 
j^— i iK(ki,k 2 ,k[,k' 2 ) 

x rrSik'jrik^k'^cs^k'jrt (2) 

where C = 7270 is the charge conjugation matrix, S is the 
fermion propagator 

S\ k a) = 772 i— a = 1,2 

iK is the interaction kernel and r a the fermion-meson 
vertex. We denote by r* its transposed and k — k u ^ u . 
The charge conjugation matrix C appears here since we 
construct the vertex function with two fermions in the 
final state. 

We have considered the following fermion (S', m) - me- 
son (0, fi) interaction Lagrangians: 

(i) Scalar coupling 

d&=gn* (3) 

for which _T Q = ig 

(ii) Pscudoscalar coupling 

dg$=ig*K><t>* (4) 

for which T a = —375. 

(ii) Massless vector exchange 

C<&=g#'fV»iP (5) 

with r a = igj^ and II ^ = —ig^v/q 2 as vector propaga- 
tor. 

Each interaction vertex has been regularized with a 
vertex form factor F(k — ft') by the replacement 

g -> gF(k - ft') 

and we have chosen F in the form: 

The BS amplitude >P is defined in terms of the vertex 
function r by: 

k 2 , P ) = s(k 1 )r(k 1 ,k 2 ,p)S(-k 2 ) 

Notice that we work with " unamputated" BS amplitude, 
which includes the external propagators S(k±), S(—k 2 ). 

Let us first consider the case of the scalar coupling and 
the corresponding ladder kernel 



(ft - ft') 2 - /i 2 + ie 
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The BS equation for the amplitude <P reads: 
i{m+\p + k) f d 4 k' 



{\p + k) 2 - m 2 + it J (27r) 



$(k',p) (7) 



{-ig 2 )F 2 (k-k') i{m-\p + k) 



(k — k') 2 — [i 2 + ie (ip — k) 2 — m 2 + ie' 

where p = ki + k 2 , k = (ki - k 2 )/2, k' = (k[ - k' 2 )/2. 

In the case of J 71 " = + state, the BS amplitude has 
the following general form: 



<!>{k,p) = Si4>i + S 2 4>2 + S 3 4> 3 + Sa4>a 



(8) 



where Si are independent spin structures (4x4 matrices) 
and 4>i are scalar functions of k 2 and p-k. 

The choice of Si is to some extent arbitrary. To benefit 
from useful orthogonality properties we have taken 

51 = 75 

5 2 = ^p 7 5 

S ' 3 = M^ 75 ~ M 75 ' 
•54 = jH<T» v p»K 75, 

where = 1(7^7^ — 7^7^)- The antisymmetry of the 
amplitude (JSJ) with respect to the permutation lf)2 im- 
plies for the scalar functions: 



01,2,4(^7^) = (/>l,2,i(-k,p) 

fo(k,p) = -fo(-k,p) 



(9) 



A decomposition similar to ((5]) was used in [5] to solve 
the BS equation for a quark-antiquark system but the so- 
lution was approximated keeping only the first term S\(j)i. 

Substituting ([5]) in cq. ([7]), multiplying it by Si and 
taking traces we get the following system of equations for 
the invariant functions <f>i. 

^ k ' P "> = [( |p + fc)2 _ m 2 + ie ] [( l p _ fc)2 _ m 2 + j £ ] 



d 4 fe' (-ig 2 )F 2 (k-k') 



y] Cgg' (fc, fc', p)(ft tt / (fc', P) (10) 



(2tt) 4 (fc-fc') 2 -M 2 + ie ^ 

a' — 1 

For scalar and pseudoscalar couplings, the coefficients 



are given by: 



Ni 



-Tr 



Si (i + k + rr?j r 1 s' J r 2 f 



k 



(11) 



Tr[S 2 } are the normalization factors, S'j 



where 2V, 

is obtained from Sj by the replacement k 



k . r r \ 



U c r^U c = r a with r a given after eqs. © and (g]). 

For the fermion-fermon vector massless coupling, the 
Cij coefficients are also given by eq. (fTTj) with the replace- 
ment J 1 ! — > gjv and F 2 — > —57" which must be con- 
tracted. For the positronium case (fermion-antifermion) , 
they are multiplied by an extra factor — 1. 



The explicit expression of for a J 71 
the scalar coupling ([3]) are 



+ state and 



1 



Cll 

C12 

Cl3 
C14 

C21 



-A'/ 



4 

mM 
C31 = 
-b'M 2 
mM 

m 2 + \ 
4 

c 2 3 = -2&'(p-fc) 
c 24 = -2b'mM 
c 32 = 2(p-fc) 
6' 



C22 



M 2 



c 33 = j [ m 



-M 
4 



2{k-pf 
M 2 



>-fc) 2 
' Af 2 



C34 



^0 m , , . 



c 4 i = M 
C42 = 2mA/ 
b' m 
C43 = 2 lM M 



C44 



where 



T 

6= a7^ 



-A/ 



M 2 fc 2 



(p-fc) 2 

(p-k)(p-k') - M 2 (k-k') 



The coefficients c^- for the pseudoscalar exchange (HJ) 
are simply obtained by changing the sign of the j ' = 1 and 
j = 4 matrix elements in the scalar ones: 



r-PS 



-C~_ 



ij=l,4 



(12) 



For the massless vector exchange ([SJ the coefficients 
are given by 



( A 


-2 








4 


-2 


-2 








-2 


-2 





U 


-2 


-2 






(13) 



The equation ([TU]) is the Minkowski space BS equation 
for two fermions we aim to solve. As in the scalar case, 
the BS amplitudes fa are singular and a direct solution of 
(fTUl) is not suitable even for the simplest kernels. 

To overcome this difficulty, the first idea is to repre- 
sent each of the BS components <pi(k,p) by means of the 
Nakanishi integral 



Uk,P) = dz' d 7 ' 

J-l Jo 



(14) 



<?i(7',*') 



[k 2 + p ■ k z' + \M 2 - m 2 - 7' + ie] 
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and apply the light-front projection to the set of coupled 
equations for the corresponding weight functions (7,(7,2). 
As mentioned in the Introduction, this projection, which 
is an essential ingredient of our previous works [TirPI] , con- 



sists in replacing k — > k 



in eq. (|10[) and integrating 



over /? in all the real domain. This integration is quite 
similar to the case of the spinless particles explained in 
detail in [I]. 

We obtain in this way a set of coupled two-dimensional 
integral equations which can be written in the general 
form valid for all types of couplings and states: 



dj'- 



1 



/•OO f+1 

Y,j o d 7 'y i dz'^( 7 ,z, 7 ',z')5 J (7'^') (15) 



where 



and 



D (7, z) = 7 + m 2 z 2 + (1 - z' z )k 



J2\ ,2 



2 2 

K = m 



M 2 



The left hand side of (|T5|) is the same for all the amplitudes 
and coincides with the scalar case pQ. 

The interaction kernel Vij can be written in the form 

V-li z V z') - { Wij (7 ' Z ' 7 '' ^ if Z ' ~ Z 

V n [J, Z, 7 , Z ) - < ^ w ^ _^ _ z ,- ) i{z > >z 



with 




(16) 
(17) 



and 



w i3 (n,z,i,z') 



am 2 (1 — z) 
~2n D 



Za{i, z, 7', z'\ v) Cij (7, z; v) j^dv (18) 



The matrix coefficients dj couple the different spin com- 
ponents. They are given in appendix [X] for the = + 
states. Notice that setting C,j = S.^ and Za = 1 the set 
of eqs. (fT5"|) decouples and each of them is identical to the 
spinless case one given in pQ . 

It turns out that, in contrast to the scalar case, most of 
the kernel matrix elements are discontinuous at z' = z. 
In some cases - like e.g. V14 displayed in ng.Q]- the value 
of the discontinuity, although being finite at fixed value 
of z, diverges when z — > ±1. This creates some numer- 
ical difficulties when computing the solutions ^(7,2:) in 
the vicinity of z = ±1. They can be in principle solved 
by properly taking into account the particular type of di- 
vergence. However the latter may depend on the particu- 
lar matrix element, on the type of coupling, the quantum 
number of the state and other details of the calculation. 
We propose in next section an alternative, and efficient 
way to overcome this difficulty. 




H=0.15 A=2 B=0.2 L=1.1 

7=0.75 Y=0.50 



Fig. 1. The kernel matrix element Vi4( 7 , z; 7', 2'), eq. (|16p . for 
z = 0.7,0.9,0.95 as a function of z' and fixed values of 7,7'. 
The discontinuity is finite at a fixed value of z but diverges 
when z 1. 



We use here a = whereas the denominator 

4-7T 1 



+ vm 2 [(1 - v)(l - z')z 2 + vz' 2 {\ - z) 



D(7, z, 7', z', v) = v(l - v)(l - z')7 
n 2 

+ vk 2 (1 - z)(l -z')[l + z- v(z - z') } 
+ (1- z)[{l-v)pi 2 +vj'} 

is the same than in our previous work [T]. 
The function 

Z A = (l-v) 2 (l-z) 2 

[3D + (1 - v)(l - z){A 2 - li 2 )] {A 2 - fi 2 ) 2 
X [D + (1 - «)(1 - z){A 2 - fi 2 )} 3 

contains the dependence on the vertex form factor param- 
eter A with 

lim Za = 1 



3 Regularized kernels 

To deal with more regular kernels, the original BS equa- 
tion (fTU|) is first multiplied on both sides of by the factor 



(m 2 - L 2 ) 



(m 2 - L 2 ) 



[(£ + k) 2 -L 2 + ie] [(f - fc) 2 - L 2 + it 



(19) 



This factor has the form of a product of two scalar prop- 
agators with mass L. It plays the role of form factor sup- 
pressing the high off-mass shell values of the constituent 
four-momenta k 2 2 = (§±fc) 2 and tends to 1 when L — > 00. 

The light-front projection and Nakanishi transform are 
then applied to the equation 



r)(k,p) (j>i{k,p) 



(20) 
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y(k,p) 

[(§ + k) 2 -m 2 + ie] [(§ - k) 2 -m 2 + ie] 

x fdV ig 2 F\k-k') A 

Since r](k,p) ^ 0, the equation thus obtained is strictly 
equivalent to (ITU1) . We will see however that, due to the 
presence of the rj factor, the light front projection modifies 
the resulting kernels which become less singular functions. 



and 

H = 1 i -^-^+ 1 ' + (l-z')(l + z) K 2 

+ (z'-z(l-z'))m 2 + ^-^-L 2 
\ / (1 z) 

To avoid spurious singularities in (|23[) due to the r] 

factor (|19f) . L 2 must be larger than to 2 , what is 

fulfilled for L > to. In practical calculations we have taken 
to = 1 and L = 1.1. 



H=0.15 B=0.2 L=1.1 
7=0.75 Y=0.50 




L n=0.15A=2B=0.2L=1.1 
7=0.75 y=0.50 



Fig. 2. The left hand side kernel V 9 (7 = 0.75, 2 = 0.7,7 
0.5, z'), eq. ()22p . as a function of z' for L=l.l m 



Fig. 3. Regularized right hand side kernel 1/14(7, z; 7', z'), ec l- 
pi) , v.s. fixed values of 2 = 0.7, 0.9, 0.95 and L = 1.1m. 
Compare with fig. [1] 



The 1 technical details of the light-front projection arc 
similar to those given in ref. [T|. The differences, due to 
the factor r]{k,p), are explained in Appendix [Bj The new 
set of equations has the following form: 



d 7 ' / dz'V 9 ( 1 ,z; 1 ',z')g t ( 1 / ,z') = 



E 



d 7 '/ d2'^( 7 ,*; 7 V) ft ( 7 V) (21) 



The kernels V 9 and V$ depend now on the parameter 
L. Closer is L to to, smoother is the kernel and more stable 
arc the numerical solutions. However the weight functions 
<7i(7, z) as well as binding energies provided by (f2~Tj) are 
independent of L. 

Notice that, in contrast to (fT5)) . the left hand side of eq. 
(|21[) is also a two-dimensional integral. The corresponding 
kernel V 9 is the same for all the components and has the 
form: 

V \n,z,1,» ) ~ p^-^^-z') if z'>z [ZZ) 



with 



W 9 ( 7 ,z; 7 ',z') 



(L 2 - TO 2 ) 
IP 



(23) 



The kernel V g is finite and vanishes for z = ±1. For a 
fixed values of 7, z and 7', is a continuous function of z' 
with a discontinuous derivative at z' = z. It is represented 
in fig. OH for z = 0.7. 

The right hand side kernel is given by: 

if z' < z 



V d Hz-,'z') = [ W fa>*->T'>* n > 
W$ is represented as a sum of two terms: 



Wt 



W a - + W b - 

13 13 



which can be written in the same form than (Tl8 



For the first term (y = a) one has 



(1-z) [D + (L 2 -m 2 )} 
L 2D Q + (1 - z)(L 2 - m 2 ) 
Z\ = Za 

Ca s~i 
13 ~ °ij 

D n = D 



(24) 



x[ Z» A ( 7 ,z, 7 ',z';v)CZ(j,z;v)^dv (25) 
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where Za , Cij and D are the same than in eq. (|T5)l . 
For the second term (y = b) one has 

Fb _ A>a - *) 

L 2D + (1 + z)(L 2 - to 2 ) 
Z\ = {l-v) 2 {l-z) 2 

[3D b + (1 - v)(l - z)(A 2 - fj, 2 )] (A 2 - /i 2 ) 2 
X [D b + (l-v)(l-z)(A 2 -fi 2 )]3 

D b = D a + v(l - v)(z - z')(L 2 - m 2 ) 

with the coefficients C\j given in Appendix [Cl 

Notice that, like Za, Z a — > 1 when A — > oo. Notice 
also that F£ — > 1 and F| — > in the limit L — > oo. One 
has consequently W a = W, W b = and V d = V given 
in (|16[) . Furthermore the peak at z' = z in the left kernel 
V 9 displayed in fig. [3] tends to a delta function S(z' — z) 
with a coefficient which reproduces the left-hand side of 
cq. (|T5|). 

For a finite value of L, both systems of equations are 
also strictly equivalent to each other but the z'-dependence 
of the regularized kernels is much more smooth and there- 
fore better adapted for obtaining accurate numerical so- 
lutions. In fig. [3] we plotted the regularized kernel VjJ as a 
function of z 1 for the same arguments 7, z, 7' and param- 
eters than in fig. [TJ where it was calculated without the 
r](k,p) factor. As one can see, the kernel is now a contin- 
uous function of z' . A discontinuity of the first derivative 
however remains at z' = z. 

Some additional remarks concerning the regularization 
procedure presented in this section are in order: 

(i) The most singular kernel is V23 . It is the only matrix 
element that after regularization by means of the rj(k,p) 
factor remains discontinuous at z' = z . The degree of the 
discontinuity in the limit z — > ±1 is however reduced and 
does no longer spoil the numerical accuracy. 

(ii) This procedure can be improved by choosing for 
rj(k,p) a more strongly decreasing function of the argu- 
ments (I ± k) . The simple replacement rj(L) — > r\(L\)r\{L<i 
would reduce the degree of singularity of the remaining 
singular matrix elements in the limit z — > ±1. We have 
not found any need for that in the present calculations. 

(in) We would like to emphasize again that despite 
the fact that the non-regularized and regularized kernels 
are very different from each other (compare e.g. the figs. 
Q] and [3]) and that the regularized ones strongly depends 
on the value of L, they provide - up to numerical errors 
- the same binding energies and weight functions ^(7, z). 
We construct in this way a family of equivalent kernels. 



4 Numerical results 

The solutions of cq. (f?T|) have been obtained using the 
same techniques than in ref [T], i. e. spline expansion of 
the Nakanishi weight functions 

9%{l,z) = ^2g i , a pS a ('Y)Sp(z) 



Table 1. Coupling constant g 2 as a function of binding en- 
ergy B for the J = state with scalar (S), pseudoscalar (PS) 
and massless vector (positronium) exchange kernels. The ver- 
tex form factor is yl = 2 and the parameter of the r\ factor 
L = 1.1. 





S 


PS 


positronium 




0.15 


0.50 


0.15 


0.50 


0.0 


B 


9 2 


9 1 


g 1 


9 2 


9 2 


0.01 


7.813 


25.23 


224.8 


422.3 


3.265 


0.02 


10.05 


29.49 


232.9 


430.1 


4.910 


0.03 


11.95 


33.01 


238.5 


435.8 


6.263 


0.04 


13.69 


36.19 


243.1 


440.4 


7.457 


0.05 


15.35 


39.19 


247.0 


444.3 


8.548 


0.10 


23.12 


52.82 


262.1 


459.9 


13.15 


0.20 


38.32 


78.25 


282.9 


480.7 


20.43 


0.30 


54.20 


103.8 


298.6 


497.4 


26.50 


0.40 


71.07 


130.7 


311.8 


515.2 


31.84 


0.50 


86.95 


157.4 


323.1 


525.9 


36.62 



on a compact domain Q = [0,7 ma J x [— 1, +1] and vali- 
dation of the equations in some well chosen points in the 
iV 7 x N z corresponding intervals. The unknown coefficients 
gi y &p and the total mass M are obtained by solving a gen- 
eralized eigenvalue problem 

XN(M)g = A(M)g 

As in the scalar case, the discrctized left hand side kernel 
N has very small eigenvalues which make difficult the so- 
lution using standard methods. This is avoided by adding 
a small constant term of the form 

N = N + e R I 

where I is the identity operator in the spline basis. The 
error in the eigenvalues thus induced is of the order of en 
and we have taken e# = 10~ 6 . 

We have computed the binding energies, defined as 
B = 2m — M, and BS amplitudes, for the J = + two 
) fermion system interacting with massive scalar (S) and 
pseudoscalar (PS) exchange kernels and for the fermion- 
antifermion system interacting with massless vector ex- 
change in Feynman gauge. In the limit of an infinite ver- 
tex form factor parameter A — > 00, the later case would 
correspond to positronium with an arbitrary value of the 
coupling constant. All the results presented in this section 
are given in the constituent mass units (m = 1) and with 
L = 1.1. 

The binding energies obtained with the form factor 
parameter A = 2 are given in table [TJ For the scalar and 
pseudoscalar cases, we present the results for fi = 0.15 
and n = 0.50 boson masses. They have been compared 
to those obtained in a previous calculation in Euclidean 
space [TBIHB] using a slightly different form factor. Notice, 
that contrary to what is written in eq. (16) of j!5) - which 
coincides with our form factor (|6]) - the vertex form factor 
used in these calculations is [12] 

A 2 

A 1 — q l + le 
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2 

g 

Fig. 4. Binding energy for scalar exchange v.s. g 2 for A = 2, 
L = 1.1, n = 0.15 and /x = 0.5. 

Once taken into account this correction, our scalar results 
are in full agreement (four digits) with [TMT6] . The pseu- 
doscalar ones show small discrepancies (s=s 0.5%). We have 
also computed |17j the binding energies by directly solv- 
ing the fermion BS equation the Euclidean space using 
a method independent of the one used in (T5]. Our Eu- 
clidean results are in full agreement with those given in 
the tabled] 




2 

g 



Fig. 5. Binding energy for pseudoscalar exchange v.s. g 2 for 
A — 2, L — 1.1, fj, = 0.15 and fj, = 0.5. 



The B(g 2 ) dependence for the scalar and pseudoscalar 
couplings is plotted in figs. @] and [5j Notice the different 
g 2 scales of both dependences. The pseudoscalar binding 
energies are fast increasing functions of g 2 and thus more 
sensitive to the accuracy of numerical methods. This sharp 
behaviour was also exhibit when solving the corresponding 
Light- Front equation (T5] . 

It is worth noticing that the stability properties of the 
BS J 7 * = + solutions for the scalar coupling are very sim- 



ilar to the Light-Front ones. In the latter case, we have 
shown jT9ll20] the existence of a critical coupling constant 
g c below which the system is stable without vertex form 
factor while for g > g c , the system "collapses", i.e. the 
spectrum is unbounded from below. The numerical value 
was found to be a c = 3.72 [50], which corresponds to 
g c = 6.84. Performing the same analysis than in our pre- 
vious work - eq. (71) from [2U] - we found that for BS 
equation the critical coupling constant is g c = 2ir, in agree- 
ment with |15j . The 10% difference between the numerical 
values of g c is apparently due to the different contents of 
the intermediate states in the two approaches. The ladder 
BS equation incorporates effectively the so-called stretch- 
boxes diagrams which are not taken into account in the 
ladder LF results. 

The positronium case deserves some comments. First 
we would like to notice that in our formalism, the singular- 
ity of the Coulomb-like kernels in terms of the momentum 
transfer l/(k — k') 2 is absent. This is a combined conse- 
quence of the Nakanishi transform ([T]) - which allows to 
integrate over k' analytically in the right hand side of the 
BS equation (|T0| - and of the consecutive light front pro- 
jection integral. After this integration, the Coulomb singu- 
larity does not anymore manifest itself in the kernel. This 
can be explicitly seen in the kernel of the Wick-Cutkosky 
model obtained in eq. (22) of our previous work [T]. 

A second remark concerns the A dependence of the 
positronium results. Using the methods developed in [191 
|2"0] we found that in the BS approach with ladder kernel 
there also exists a critical value of the coupling constant 
g c = ir. Note that, as in the scalar coupling, the very 
existence and the value of this critical coupling constant 
is independent on the constituent (m) and exchange \x 
masses but depends on the quantum number of the state. 

Table 2. Coupling constant g 2 as a function of binding energy 
B for the positronium J = state in BS equation in the region 
of stability without vertex form factor (A — > oo), i.e. g < n. 
They are compared to the non relativistic results. 



B 


2 

9nr 


„2 

9bs 


0.01 


2.51 


3.18 


0.02 


3.55 


4.65 


0.03 


4.35 


5.75 


0.04 


5.03 


6.64 


0.05 


5.62 


7.38 


0.06 


7.95 


8.02 


0.07 


11.24 


8.57 


0.08 


13.77 


9.06 


0.09 


15.90 


9.49 



The ground state positronium binding energies with- 
out vertex form factor are given in table [5] for values of 
the coupling below g Cl Nonrelativistic results 

9nr = Stt^B/ui 

are included for comparison. One can see that the rela- 
tivistic effects are repulsive. 
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These results arc displayed in fig. [5] (black solid line), 
and compared to the binding energies obtained with two 
values of the form factor parameter A = 2 (dashed) and 
A = 5 (dot-dashed). The stability region is limited by a 
vertical dotted line at g = g c = n. Beyond this value the 
binding energy without form factor becomes infinite and 
we have found B(g — > g c ) w 0.10. 



B/m o.i 




Fig. 6. Binding energy for J = positronium state versus g 2 
(black solid line) in the stability region g < g c — n. Dashed and 
dotted-dashed curves correspond to the results for increasing 
values of the vertex form factor parameter A = 2 and A = 5 
respectively. They are compared to the non relativistic results 
(dotted line). 



The inclusion of the form form factor has a repulsive 
effect, i.e. for a fixed value of the coupling constant it pro- 
vides a binding energy of the system which is smaller than 
in the A —> oo limit (no cut-off). This is also illustrated in 
fig. [7] where we have plotted the A dependence of g 2 for 
two different energies. One can see that the value of the 
coupling constant to produce a bound state is a decreasing 
function of A. The size of the effect depends strongly on 
the binding energy but for both energies the asymptotics 
is reached at A w 20. This behaviour is understandable in 
terms of regularizing the short range singularity of — 1/r 2 
interactions. 

Finally we present some examples of the Nakanishi 
weigh functions (^(7, z). They correspond to a B = 0.1 
state with the scalar coupling and the same parameters 
A = 2, fj, = 0.50 than in table [TJ They are displayed in 
fig. [5J In the upper figure is shown the 7-dependcnce for a 
fixed value of z and in the lower figure the z-dependence 
for a fixed 7. Notice the regular behaviour of these func- 
tions as well as their well defined parity with respect to z 
~ 91,2,4 are even and g% is odd - consequence of relations 
© . As in the scalar case the e^-dependence of gi is more 
important than for the binding energy. 

Corresponding BS amplitudes 4>i are displayed in fig. 
[5J Their computation by directly applying equation (| 14|) 
is not very suitable due to the zeroes of the integrand 





Positronium J=0 




B=0.09 1 




B=0.05 = 





5 10 15 20 25 30 35 40 45 50 

A 

Fig. 7. /l-dependence of g 2 for J = positronium state for 
fixed values of binding energies B = 0.05 and B = 0.09. 




. T&0.54 
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A=2u=0.50 B=0.10 


^*V* ■ 












E=10^ 








Fig. 8. Nakanishi weight functions for scalar exchange for A = 
2, B = 0.10, fx = 0.50 and L = 1.1. 



denominator. We have preferred to compute <f)i by first 
factorizing the pole singularities in the right hand side of 
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cq. una 

</>a{k,p) 



i 2 r a (k, P ) 



[(hp + fc) 2 — m 2 + ie][(ip — fc) 2 — m 2 + ie] 



and computing r a which is a regular function of fc. By 
inserting the Nakanishi transform ((T4)) in the expression 
for r a , the integral over k! in (| 10[) can be performed ana- 
lytically and obtains the form 

f a {k,p) = ^2j dj J dz I aa ,(k,p,"/,z)g a i(j,z) 



where I aa : is a more appropriate function for computation 
purposes. 

The upper figure [S] represents the fco dependence oi 
<j>i for a fixed value of | fc |= 0.2. They exhibit a singu- 
lar behaviour which corresponds to the pole of free prop- 
agators fco = e k — ^f- The lower figure represents the 
| fc | dependence of the amplitudes <f>i for a fixed value 
fco = 0.04. For this choice of arguments, the amplitudes 
are smooth functions of | fc |, though they will be also 
singular for fco > y = 0.05. Notice that all the infini- 
ties in the BS amplitudes <fn come from the free propa- 
gator. However t a , although being smooth, has its own 
non-analytical points. For instance it obtains an imagi- 
nary part for ± fc) 2 > (m + /i) 2 . 



5 Conclusions 

We present a new method for obtaining the solutions oi 
the Bethe-Salpeter equation in Minkowski space for the 
two-fcrmion system. It is based on a Nakanishi integral 
representation of the amplitude and Light-Front projec- 
tion and constitutes a natural extension our previous work 
for the scalar case PQ. 

A straightforward generalization of this approach, how- 
ever results into a singular fermionic kernel. In order to 
smooth this singularities, a proper regularization of the 
kernels has been proposed. This generates a family of 
strictly equivalent equations depending on one parame- 
ter L. Their solution gives the same binding energies and 
Nakanishi weight functions ft (7, z) . 

The binding energies for the scalar and pseudoscalar 
exchange kernels and for massless vector exchange (positro- 
nium) are found. They coincide with the ones found via 
Euclidean space solution, thus providing a validity test of 
our method. The solutions for the scalar and positronium 
states without vertex form factor (A — > 00) arc found to 
be stable below some critical value g c of the coupling con- 
stant, respectively g c = 2?t (scalar) and g c = tt (positron- 
ium). 

The BS amplitudes are obtained in terms of the com- 
puted Nakanishi weight functions. The exhibit a singular 
behaviour due, on one hand, to the poles of the free propa- 
gators and, on the other hand, to non analytic branch cuts 
which are responsible for the appearance of an imaginary 
part. 



L m=1 u=0.50 B=0.1 / / 


1 1 E 


L k=0.2 / / 




j 


J 






L °; 




a , 1 , 1 , 1 





m=1 [i=0.50 B=0.1 

k =0.04 




Fig. 9. Bethe-Salpeter Minkowski amplitudes corresponding 
to fig. [8] Here fc =| fc |. The amplitudes <j>\ and 4>2 are indis- 
tinguishable. 



An alternative approach to avoid the singularities of 
the Minkowski amplitude would be to solve the BS equa- 
tion in the Euclidean space in terms of the Nakanishi 
weight functions. Once solved, we can again restore the 
Minkowski BS amplitude by the integral ([T]). These differ- 
ent methods should provide the same solutions although 
with different precision and stability. The results of this 
study will be published in a future paper |17j . 
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A The coefficients CV, 

The spin-coupling matrix elements (7, z; v) in cq. (|18l) 
determining the kernel of the equation (fT5j) . and also the 
kernel (j25)) at v = a, for the J = state and scalar cou- 
pling are: 



C12 
C 13 

C14 
C21 
C22 
C23 
C24 
C3I 
C32 

C33 
C34 
C41 
C42 
C43 

C44 



(I-2) 4m 2 

M 

4m 



(1— v) K+K- 
(1-z) 2 4m' 2 M 2 

Cl2 

(7+m 2 ) 2 , 
(l-z) 2 2m ! M 2 " t " 

l-t. g+ if- Q 



72 + (2-z)m 2 
4m 2 (l-z) 



(l-z 2 )M 2 
32m 2 



2m 2 A/ 4 
if + if _ 



= 



' (1-z) 2 2mM J 



Q 



(1-z) 2m 2 

1-u Q- Q+ 
(1-z) 2 2m 2 A/ 2 
!-■» Q- 
1-z 2mA/ 
Af 2 
4m 2 

2C 12 
C34 

l-v 47-4(l-22)m 2 + (l-z) 2 M 2 
1-z 16m 2 



with 



/ M 
K± = 7 + I m± y(l 

Q± = 7 + ™ 2 ±— (1 



-) 2 



The coefficients for the pseudoscalar and massless vec- 
tor exchange are given by the relations (|12[) and (fT3| re- 
spectively. 



B Derivation of the equation (I2T1) 

In order to derive the equation (|21[) . we substitute in eq. 
(|2T)|) the amplitudes 4>i(k,p) in the form (Q]), shift the vari- 
able k — > k + -^(3 and integrate over (3 in the infinite 
limits. 

The integrand depends on the scalar variables k 2 , p-k 
and — which are not independent but are related to each 

other by a constrain. Following appendix A from [1], we 
parametrize them in terms of two independent variables 

7, z as: 



k 2 

uj-k 
oj-p 

p-k 



(7 + z 2 m 2 ) 
1-z 2 



s[7 + z 2 m 2 + (1 — z 2 )k 2 



where k 2 



i 2 — jM 2 . In the spinless case [I], when we 
applied integration over (3 to the equation without the 



factor 7](k,p), we found that its left-hand side obtained 
the form: 



l.h.s. 



dz' I d 7 ' / d{3 

J — OO 

9(l', z') 



(26) 



f3{z - z>) + i + (1 - zz>) (k 2 + 



As a function of /3, it contains one pole at 



(27) 



Taking the residue, we found that the result of integration 
is proportional to S(z' — z) and consequently the left-hand 
side obtained the form of a one-dimensional integral: 



l.h.s. 



iirg(j',z)d^ 



„,/ 



7' + 7 + z 2 m 2 + (1 — z 2 )k 2 



The factor rj(k,p) itself consists of two factors. There- 
fore, after replacement k — > k + j3, the factor rj(k + 
/3, p) v.s. (3 obtains two poles: 



n = 



1 - z 



+ 



l + z 



Integrating then the left-hand side of (|20|) (with shifted 
argument k + -^-0) over /?, we can take into account one 
of these poles, which is in the opposite half-plane than 
the pole (|2"T)) . The position of the pole (|2"T)) depends on 
the sign of the difference z — z' . The integration, after 
cancellation of both sides of equation by a common factor, 
results into the left-hand side of eq. (f2~Tj). Now the delta- 
function 6(z' — z) does not appear and the integral is two- 
dimensional one. 

Concerning the calculation of the right-hand side, the 
difference relative to what was explained in [T] for the 
spinless particles is minimal. In the spinless case, the right- 
hand side already had two propagators, like in eq. (fTU)) . 
Their poles v.s. (3 (as well as the pole of the kernel) were 
taken into account. The factor r](k,p) adds two new poles 
in other points, which differ from the scalar case by the re- 
placement m — y L and which can be taken into account in 
a similar way. After shifting the variable k — > k+ ^ (3 the 
coefficients c,j in right-hand side of (|10|) become depend 
on (3, but they have no poles. 

In this way, multiplying both parts of eq. (|10p by r](k,p), 
integrating analytically over k'), shifting the four- vector 
k — > k + -^j- f3 and integrating over /?, we find the system 
of equations (j2"Tj) . 
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C The coefficients C, 



The spin-coupling matrix C*-(7, z; v) in the kernel (|25[) at 
v = b for the J = state and scalar coupling are: 



t~<a 
°12 





Cl2 
rib 
°13 

(~*b 

ls u — 

Cb r^a 

21 ~ °21 

(~*b r^a 

°22 — °22 

r^b 

°23 — 

C b 24 = 

C 3 \ = 



r^b 
°32 

C*b 

°33 

(~*b 
°34 

°41 

/~*b 

°42 

°43 
°44 



33 



°41 

2C 12 

r*b 
°34 

°44 



z(L 2 -rn 2 ) 



(1-z) 2 4m 2 Af 2 

(L 2 -m 2 ) [47+2m 2 +2L 2 ~(l-z)Af 2 ] 
4(l-z) 2 m 2 A/ 2 
1-v K b + K b _Q"_ 
(1-z) 3 2m 2 M 4 

(1-z) 2 2mA/ 3 



Q- 

z) 2m' 2 



(1-z) 

1 



-m 2 ) [47+2m 2 +2L 2 -(l-z)A/ 2 ] 



4m 2 M 2 



2mM 



°42 



(1— i>)z(L 2 -to 2 ) 
4m 2 (l-z) 



12. 
13. 

14. 

1."). 

16. 
17. 
18. 

10. 

20. 



N. Nakanishi, Phys. Rev. 130, 1230 (1963); Prog. Theor. 
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with 



Q-=i 



2 1 



9 M 2 , 



and C§ 



Cij given in Appendix |A"1 



The relations with the pseudoscalar and massless vec- 
tor exchange are also given by (fT2"j) and (fT5)) respectively. 
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